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IK AIRPLANE DYNAMICS 
By Robert T. Jones 

SUMMARY 



The method of operators is used in the application of 
nonunif orm-lif t theory to problems of airplane dynamics. 
The method is adapted to the determination of the lift under 
prescribed conditions of motion or to the determination of 
the motions with prescribed disturbing forces. 



INTRODUCTION 



Problems in airplane dynamics are usually treated on 
the assumption that the air forces are instantly adjusted 
to each motion of the airplane. Since the development of 
recent theories for the nonuniform motion of airfoils, it 
has become possible to consider more exact laws f or the ad- 
justment of the lift. 

The nonunif orm-lif t theory has already been applied to 
certain dynamical problems, notably to the problem of flut- 
ter. These applications have, however, been confined either 
to approximate solutions or to cases in which the t y p e o f 
motion is prescribed beforehand. The more usual problem, 
in which the resulting motion is unknown, requires tho so- 
lution of integral equations. The present paper shows ho_w 
solutions of these equations may be obtained fairly simply 
by operational methods. 



SUPERPOSITION OP LIFTS 



In nearly every aerodynamic problem, the approximations 
that must be made to effect solutions are such as to lead to 
linear relations. Thus, in the case of the unsteady lift of 
a wing, Laplace's equation combined with the assumption of 



2 



N.A.C.A. Technical Note No. 667 



an undistorted wake leads to a linear relation between the 
lift and the angle of attack. Such a relation means that 
the lift due to the sum of" two variable motions is equal 
to the gum of the lifts for the two motions taken independ- 
ently. 

In particular, if the lift following a sudden unit 
jump of angle of attack, is known (see .reference 1), then 
the lift for any variable motion is easily obtained by 
breaking the given motion down into a succession of small .. 
jumps or steps and adding the lifts incident to each one. 
The case treated by Wagner thus becomes the key to the 
calculation of lift for any variable motion. 



Wagner's function (reference 1) giving the lift after 
a sudden unit jump of angle of attack ( t wo" - dimens iona l 
case) may be denoted by c, (s). The superposition of 

lifts for any variable motion a(s), as previously ex- 
plained, is accomplished by the integration of JDuhamol's 
integral 

s 

/ c (s - s 0 ) a 1 (s 0 ) ds 0 (1) 

'o li ~ ~ : 

(See reference 2.) 



OPERATIONAL SOLUTION OF INTEGRAL EQUATIONS 



It is evident that, in order to take account of— un- 
steady air-flow phenomena in the theory of airplane dynam- 
ics (including stability and related problems) the custom- 
ary instantaneous equations of notion must be replaced by ■■ 
equations involving the integral (1), The equations of 
jnotion then become linear integral equations. Solutions 
of these equations may be conveniontly obtained by opera- 
tional mothods. 

Let D represent the operator d/ds and let 1 = l(s) 

represent the unit jump fuction, that— is, a function of s_ 

having the value 1 at s > 0 "and having the value 0 at" 

s < 0, Then a function of s may be represented by a com- 
bination of operations on the unit jump function 



cp(s) =^(D)lc>) 



(2) 
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The combination, of ' operations cp(D) on 1 necessary 
to reproduce the function cp(s) is called the "operational 
equivalent" of the function" cp(s). The operational equiva- 
lent ip of a given function cp may be found by the infi- 
nite-integral theorem (reference 2) 



cp(a) = a J cp(x) e~ ax dx (3) 
o 

A general operational equivalent is 

s n = T (1 + n) D" n l(s) (4) 

(3oe Peirce's table, p. 63, HO. 493.) 

Tho operational treatment of integral equations is 
based on the proposition that an integral of the form 



<p(s) = Z(s) X(O) +J* Z(s - s 0 ) X'(s 0 )-ds 0 (5) 

maj' be regarded as the solution of a linear differential 
equation. As such, its operational equivalent is 

(p (s) = Z(D) X(s) = Z(D) X(D) 1 (s) (6) 



wiiere Z and X are the operational equivalents of the 
functions Z and X. 

In order to illustrate the operational solution, let 
it be required to find the function X(s) from 

s 

/ , ds 0 - 2X(s) = — s - s (7) 

^o - s 0 



4 



N.A.C.A. Technical Koto No; 667 



assuming that Z(0) = 0. Here the function corresponding 
to Z(s) in equation (5) is lA/s. With the aid of equa- 
tion (4), the various componenet.s are written in operational 
form ■ , ... 



Z(s) = rd/2) 1 (e) (8) 



Equation (6) "becomes 

T(l/2) X(s) - 2X(s) = ~ L) 1 ( s ) 

r (i/2) 2 

or X(s) = l(s) 

ill/2) V^D - 2 • 



Simplifying: 2 

r(i/2) 

1 ^~ 1 1 
x( s ) = i = 4r* = * s s (io) 



D 



a 2D 2 



r(i/2) -■ — 



It will be helpful to review, certain aspects of the 
theory of unsteady lift before proceeding to- -the applica- 
tion of this theory in airplane dynamics* It is found 
convenient to think of the lift on the airf-oirl as composed 
of three parts; (1) A part due to instantaneous acceler- 
ation of the noncirculatory potential flow. This lift is 
equal to the virtual additional mass of the wing 

c a ' - - " ~ ' 

(tt — P per unit span, for infinite aspect ratio"* times 

the rate of increase of the relative wind velocity normal 
to the chord. (2) A part due to the .circulatory flow 
and dependent on the angle-of -at tack variation, i.e., the 
lift ^iven by c^ (s). (S) A part due to the circulatory 

l 
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flow and ascribed to a relative curvature or camber of 
the airfoil in pitching motion. The third component will 
"be automatically included with the second if the angle of 
attack is obtained "by resolving velocities at the 75-per- 
cent-chord point. 

With these provisions, the instantaneous lift of an 
airfoil in combined pitching and vertical motion may be 
written (see fig. l) 

o l ( 8 ) = p, fa' (s) + i Bo &" (s) J 

_ O£f0~* , * ^ 

+ [a(0) + l 7B 9» (0) o l (s) + / c L (b- s 0 )U« _(b 0 ) + l 7S 6 ,, (s 0 5 j ds 

°" L (11) 

where [i is a coefficient fof'the virtual additional mass 
of the wing ((j, = tt for infinite aspect ratio). Now let 
c^ be the operational equivalent of Wagner' s functio_n c^ 

c l ( 8 ) = UL |^Da+l eo D2eJ + c li (D) [alSp l ?5 D6fej (12) 

No concise formula for ~cj (s) is known although it 
is found that Wagner's curve is 1 reproduced almost exactly 
by the equation 

c, ( s ) = C 0 + C^e 5 ^ 3 4- C s e X2S (13) 

l _ 

where C Q = ,2tt 

G 1 = -0.330- tt 

0 2 = -0.670 tt 

X x = -0.0455 

A a =• -0.300 



and where s refers to the half-chord as unit, that is 
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Vt 
c/3 

In this form, the operational equivalent is readily 
found from the relation 

e* 8 = ' l(s) (See reference 2) (14) 

D - A 

whence 

D D 

cj (D) = C 0 H- C + 0 3 (15) 

1 D - Xl D - Xa 

The calculation of lift under a prescribed variation 
of angle of attack can "oe illustrated— by assuming that the 
airfoil is given a sinusoidal motion 

a(s) (or 6(e)) = E.p. or I. P. of e ins (16) 

This variation is reduced to operational form (see equation 
(14)): 

a(D) = -- - -- - (17) 

D - in 

c z ( s ) = p. D — 2 i + (o Q + 0l ' JL-+ C a _L-)_£~l< (18) 

n D - in v D-A, D-Aa^D-in 



(See equations (12) and (15),) 

The resulting operator may be- evaluated by the Hoavi- 
side expansion theorem: 



f(3) _ f (0) 2 f (X) As 
F(0) 1 ~ F(O) + A A3P 1 (A) 6 (19) 

where the A's are the roots of S'(-D) = 0. 
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o <s) = p. ine ills + To o + O ir ^- + C a -^-T 
l a L in - X x in - X a J 



ins 



+ 0,—-^— e AlS + C 2 — 2 e*s s . (20) 

Ai - in X 2 - -in 

X s 

The terms involving e disappear in time and henco 
may "be disregarded in. a continuous oscilla.tion. The terms 

e ine (21) 



yield approximate expressions of the lift functions for 
tno oscillating airfoil introduced by Thoodorsen fref er- 
ence 3) : "" " 

2 a 
n n 




2 T TS , (n) = C Q + O l + Cg 

A a + n 

(22) 



, a s ■ . s a 

A x + n X 2 + n 



X x n A 2 n . 

2ttCt (n) = - C,— — - C a — 7 + r v\ - ~ " 

A x + n A 2 + n / 

As pointed out by Garrick (reference 4), Theodrosen's 
function for sinusoidal motion 



C(in) = F(n) + iG(n) (equation (21)) 

may bo regarded as the operational equivalent of Wagner's 
curve , i.e., 

2ttO(D) 1 = c l (D) 1 = G l (s) (23) 



This fact may be verified by referring to equation (15). 
This relation is especially interesting because it shows a 
connection between the Fourier and the operational analyses 
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Thus, if the response of a linear system to a continuous 
sinusoidal excitation is known, 

S n (s) = f (in) o ins (24) 



Then the function f furnishes immediately the operational 
equivalent of the unit-response so that for any variable 
exci tat i on Z ( s) , 

R(s) = f (D) Z(s) = f (D) Z(D) 1 (25) 



In general , the motion of the airfoil or airplane 
will not "be prescribed befx>-rehand but must be determined 
from dynamical equations. This typne- of problem can be il- 
lustrated simply by considering the disturbed vertical mo- 
tion of the airplane without pitching. The dynamical 
equation in this case is 

cL w 

m — - resisting force = impressed force, Z (26) 
dt 



m 



where w is the vertical velocity of the airplane and 
is the mass including the virtual additional . mass of the 
wing. Since 

in _ jL is (27) 

dt ~ c/2 ds {Z7) 

m dw = _2m__ x g p/2 c/3 x _f_ da (28) 
dt S p/2 c c/2 ds 

Making the ■ sub sti tut i on 



2 m 



g (29) 



S P/2 c 

and writing the equation in coefficient form, 

cr.Dct + (D) a = o l ( s ) (30) 
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where c is the lift coefficient of the given disturb- 
l o 

ing force. The operational solution is 

cl(b) = ~ o l (s) (31) 

crD + c ti (D) o 

Again, as in the case, of the lift, the solution for the 
elementary jump is the key to solutions for variable con- 
ditions. 

1 

CL (s) = 1 (32) 

<TD + c, (D) 

Replacing c. (D) by (15) and simplifying: 

(D - X x ) (D - X s ) f(D) , 
a x (s) = 1. = 1 (33) 

aD 3 + bD a + CD + d P(D) 

which is in standard form for evaluation by the expansion 
theroem (19). Finally, 

a(s) = c (o) a,(s) + f a (s - s 0 ) c ' (s 0 ) ds 0 (34) 
l ° Jo 0 

The extension of this treatment to problems involv-_ 
ing a number of degrees of freedom will b© evident. 



Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronaut ics, 
Langley Field, Va. , September 12, 1938. 
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Figure 1.- Moving aj-es. a = w/v ; s = Vt^c/i 



